Electronic structure of magnetically modulated graphene 
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'^ ' Abstract 

CN ' We present a theoretical study of the electronic structure of magnetically modulated graphene. 

We consider monolayer graphene in the presence of a perpendicular magnetic field and a unidirec- 
r1 , tional weak magnetic modulation. The density of states and the bandwidth of the Dirac electrons 
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in this system are determined. We have found magnetic Weiss oscillations in the bandwidth and 
the density of states. These oscillations are out of phase and larger in amplitude than the ones 
in the electrically modulated graphene. In addition, these oscillations are in phase and smaller in 
amplitude to those of magnetically modulated standard electron gas system. 



I. INTRODUCTION 

The successful preparation of monolayer graphene has allowed the possibility of study- 
ing the properties of electrons in this systeni[l|. The nature of quasiparticles called Dirac 
electrons in these two-dimensional systems is very different from those of the conventional 
two-dimensional electron gas (2DEG) systems realized in semiconductor heterostructures. 
Graphene has a honeycomb lattice of carbon atoms. The quasiparticles in graphene have a 
band structure in which electron and hole bands touch at two points in the Brillouin zone. 
At these Dirac points the quasiparticles obey the massless Dirac equation. In other words, 
they behave as massless Dirac particles leading to a linear dispersion relation ek = vk { with 
the characteristic velocity v ^ lO^m/s). This difference in the nature of the quasiparticles in 
graphene from conventional 2DEG has given rise to a host of new and unusual phenomena 
such as anomalous quantum Hall effects and a vr Berry phase [l|, |2|, y]. Besides the funda- 
mental interest in understanding the electronic properties of graphene there is also serious 
suggestions that it can serve as the building block for nanoelectronic devices j^ . Since Dirac 
electrons can not be confined by electrostatic potentials due to the Klein's paradox it was 
suggested that magnetic confinement be considered |5|. Technology for this already exists 
as the required magnetic field can be created by having ferromagnetic or superconducting 
layers beneath the substrate 161. 

In conventional 2DEG systems, electron transport in the presence of electric [7|, |8|,|9| and 



magnetic modulation [I0|,lll|,ll2] has continued to be an active area of research. In graphene, 
electrical transport, density of states, bandwidth, thermodynamic properties and collective 
excitations in t 



dictions made 



le presence of electrical modulation have been considered and theoretical pre- 
13l . 1 14 . Il5| . It is interesting to study the affects of the magnetic modulation 



on the Dirac electrons in a graphene monolayer and investigation in this direction has re- 



cently been carried out 16|, llTI] ■ In this work we study the effects of magnetic modulation on 
the bandwidth (A) and the density of states (DOS) of the Dirac electrons in graphene. Both 
these quantities are essential ingredients in studying problems such as electron transport, 
collective excitations and thermodynamic properties etc. We consider an external magnetic 
field perpendicular to the graphene monolayer that is modulated weakly and periodically 
along one direction. We find that the magnetic Weiss oscillations in the bandwidth and 
density of states are out of phase and larger in amplitude than the electric Weiss oscillations 



found in the corresponding electrically modulated system. 

In section II, we present the formulation of the problem. Section III contains the calcula- 
tion of the density of states. Bandwidth for magnetically modulated graphene is discussed 
in section IV followed by conclusions in section V. 

II. FORMULATION 

We consider two-dimensional Dirac electrons in graphene moving in the x-y-plane. The 
magnetic field {B) is applied along the z-direction perpendicular to the graphene plane. The 
perpendicular magnetic field B is modulated weakly and periodically along one direction 
such that B = (5 + Bocos{Kx))^. Here Bq is the strength of the magnetic modulation. 
In this work we consider the modulation to be weak such that Bq << B. We consider the 
graphene layer within the single electron approximation. The low energy excitations are 
described by the two-dimensional (2D) Dirac like Hamiltonian (h = c = 1 here) [l|, y, Il6l | 

i7 = i;V.(-iV + eA). (1) 

Here a = { a x, ay} are the Pauli matrices and v characterizes the electron velocity. We 
employ the Landau gauge and write the vector potential as A = (0, Bx+{Bq/K) sin(i^x), 0) 
where K = 271 /a and a is the period of the modulation. The Hamiltonian given by Eq. (1) 
can be expressed as 

H = —iv cr .V + evcTyBx + evCy — sm{Kx). 

The above Hamiltonian can be written as 

H = Ho + V'{x) (2) 

where Hq is the unmodulated Hamiltonian and V is the perturbation potential due to the 
periodic magnetic field in one direction such that 

Hq = —iv^o^.V + evo^Bx 

and 

V'{x) = ujQay sui{Kx). 



where luq = ^^^. The Landau level energy eigenvalues without modulation are given by 

e{n) = LUgy/n (3) 

where n is an integer and Ug = v\/2eB. As has been pointed out [l3|, [l6| the Landau level 
spectrum for Dirac electrons is significantly different from the spectrum for electrons in 
conventional 2DEG which is given as e{n) = uJc{n + 1/2), where ujc = eB/m is the cyclotron 
frequency. 

The eigenf unctions without modulation are given by 

^n,ky{r) = —^= I (4) 



where 



p (X Xq) / ZL ^ ^ 

$n(x,Xo) = , _ =Hn' 
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where / = \Jl/eB is the magnetic length, xq = l^ky, Ly is the y-dimension of the graphene 
layer and Hn{x) are the Hermite polynomials. Since we are considering weak modulation 
Bq << B, we can apply standard perturbation theory to determine the first order corrections 
to the unmodulated energy eigenvalues in the presence of modulation 

OO Ly 

Ae„,,^ = jdxjdy^l,^{T)H\x)^^^u^{T) (5) 



-OO 



with the result 



^^u,ky = ^0 cos(is:a;o) ( — [Ln-i{u) - L„(n)] j (6) 

where u = K'^P/2 and Ln{u) are the Laguerre polynomials. Hence the energy eigenvalues 
in the presence of modulation are 

e{n, ky) = e{n) + Ae^ ^^ = Ug^/n + iOoGn{u) cos{Kxq) (7) 

with Gn{u) = -^e~"/^[L„„i('u) — L„(m)]. We observe that the degeneracy of the Landau 
level spectrum of the unmodulated system with respect to ky is lifted in the presence of 
modulation with the explicit presence of ky in xq. The n = Landau level is different 
from the rest of the levels as the energy of this level vanishes and no modulation induced 
broadening of this level occurs. The rest of the Landau levels broaden into bands. The 



Landau bandwidths ~ Gn oscillate as a function of n since Ln{u) are oscillatory functions 
of the index n. 

Before we begin the calculation of the density of states and the bandwidth, it is necessary 
to discuss the regime of validity of the perturbation theory presented here. For large n the 
level spacing given by Eq. (3) goes as uOg{^/n — \/{n — 1)) — > ^9 27n ^"^^ ^^^ width of 
the nth level given by Eq.(6) goes as "^^^^i — • There is therefore a value of n at which the 
width becomes equal to the spacing and the perturbation theory is no longer valid. This 
occurs when n^ax = V^tt^ J- where B' = -^ = 0.0054T for a fixed value of a = 350nm. For 
a fixed electron density and the period of modulation this suggests the maximum value of 
the magnetic modulation Bq above which it is necessary to carry out a more sophisticated 
analysis. 

At this stage we can compare the energy spectrum of Dirac electrons in magnetically mod- 
ulated graphene with both the electrically modulated graphene system and the electrically 
modulated standard electron system. The differences are: 

• The standard electron unperturbed energy eigenvalues depend linearly on the magnetic 
field and the quantum number n while they depend on the square root of both the magnetic 
field and n for the Dirac electrons in graphene. 

• In magnetically modulated graphene, we have the difference of two successive Laguerre 
polynomials while for standard electrons and electrically modulated graphene, we have the 
sum and average of the two successive Laguerre polynomials respectively. 

• In magnetically modulated graphene the perturbed energy eigenvalues due to modula- 
tion are multiplied by the square root of the Landau band index ^/n that was absent in the 
expression for the electric case. 

These differences will affect the density of states and the band width as we show in the 
next section. Note that for a weak magnetic modulation case under consideration the quan- 
tum numbers n can be referred to as the magnetic Landau band indices and are equivalent 
to the Landau level quantum number n for the unmodulated system. Thus the magnetic 
modulation induced broadening of energy spectrum is non-uniform, a feature which will be 
of significance in understanding the behavior of Dirac electrons in modulated graphene. 



III. THE DENSITY OF STATES (DOS) 

It is well known that in the absence of both the magnetic field and modulation, the 
DOS consists of a series of delta functions at energies equal to e{n). The addition of a weak 
spatially periodic magnetic modulation however modifies the former delta function like DOS 
by broadening the singularities at the energies(£(n)) into bands. The density of states is 
given by 

D{e) = jY.^{e-e{n,ky)). (8) 

The sum on n extends over all occupied Landau levels and A is area of the sample. By using 
the energy eigenvalues given in equation (7), we can express D{e) as: 

a 

D{e) = 2^^ XI jd^o5{e - £„ - |G„| cos{Kxo)), (9) 

where e„ = ojg^fn, and a factor 2 is due to spin degeneracy. Evaluation of the Xq— integral 
in the above equation yields the zero temperature density of states of the density modulated 
two-dimensional Dirac electrons (DM2DDE): 

where Q[x) is the Heaviside unit step function. Here we can see that the one-dimensional 
van Hove singularities of the inverse square-root type exist on either side of the low and high 
energy edges of the broadened Landau bands forming a double peak like structure. 

The dimensionless density of state as a function of energy at i? = 0.35T is shown graph- 



th t 



le magnetically and electrically modulated 
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16[: V ^ lO^m/s, riD = 3 x 10^^ m ^, 



ically in Fig.(l) as a function of 1/-B, for bot 
graphene using the following parameters 
a = 350nm, c^o = 1 meV, and kp = (2nn£,Y^'^ being the Fermi wave number of the un- 
modulated system in the absence of magnetic field. The modulation induced pronounced 
oscillations are apparent in the weak magnetic field regime, these are the Weiss oscillations, 
superimposed on SdH-type oscillations in the high field regime that are not induced in our 
results but it is interesting to highlight their characteristics. The origin of these two types 
of oscillations can be understood by a closer analytic examination of equation (10). In the 
regime Ug > |G„|, the unit step function vanishes for all but the highest occupied Lan- 
dau band corresponding, say, to the band index N. Hence the sum over n is trivial. The 



analytic structure primarily responsible for the Weiss type of oscillations is the function 
^dCnl — \s — en\)i which jumps periodically from zero (when the Fermi level is above the 
highest occupied Landau band) to unity (when the Fermi level is contained with in the high- 
est occupied Landau band), these oscillations are with constant period in 1/B similar to 
the SdH type of oscillations in case of electrically modulated graphene but are out of phase 
and larger in amplitude. On the other hand, modulation of the amplitude of the Weiss os- 
cillations displayed in Fig. (1) is largely a consequence of the oscillatory factor Gn-, which 
exhibits commensurability oscillations. We also find that the minima of the density of states 
in the magnetically modulated system occur at the maxima of the electrically modulated 
one with the result that the oscillations in the density of states in the two systems are out 
of phase. We also observe that broadening of the Landau levels is greater in magnetically 
modulated graphene compared to electrically modulated graphene. In addition, van Hove 
singularities forming double peak like structures appear at the low and high energy edges of 
the broadened Landau levels. 

IV. THE BANDWIDTH (A) 

To better appreciate the modulation of the amplitude of Weiss oscillations we plot the 
bandwidth as a function of the magnetic field in Fig. 3. The width of the nth Landau level 
is given as 

A = 2 |G„| = l^e-"/2[L„_i(«) - LM] (H) 



This is clearly different from the electrical 



ulated standard electron result 



id, 
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ly modulated graphene and magnetically mod- 
12|. The bandwidth is plotted for n = np where 
Up = Ep/Ug is the Landau level index at the Fermi energy similar to the electrically mod- 
ulated graphene system [to compare ^p = ^ — |,with tOc = ^ for standard electrons in 
2DEG]. In Fig. (2) we present the bandwidths as a function of the magnetic field for both 
the magnetically modulated graphene and electrically modulated graphene. The parame- 
ters used in the figures are: v ~ lO^m/s, rip, = 3 x 10^^ m~^, a = 350nm, ujq = 1 meV, 
and kp = (27rn£))^/^. The bandwidths as a function of the magnetic field for the magnet- 
ically modulated graphene and magnetically modulated 2DEG are shown in Fig. (3). For 
the low magnetic fields under consideration, the magnetic modulation induced bandwidth 
in graphene is out of phase and larger in amplitude (~ 1.5 times larger at magnetic field 



B = IT) than the electric modulation induced bandwidth in graphene. In addition, it is in 



phase but with smaller amplitude (~ 5 times sma^ 
to the magnetically modulate 2DEG system. [10|, lU 



ler at magnetic field B = IT) compared 



12| . An analysis of the observed change 



in amplitude of the oscillations in the bandwidths of the three systems considered here is 
presented in the next section. 

A. Asymptotic Expression 

The asymptotic expression of bandwidth can be obtained by using the folowing asymp- 
totic expression for the Laguerre polynomials in the limit of large n as 

exp-"/^L„(M) ^ -^=^ cos(2v/^ - -) (12) 

Using n = Up = Ep/Ug and substituting the asymptotic expression given by equation (12) 
into equation (11) yields the asymptotic expression for the bandwidth 

Uo^/riF 1 



aUjQ^Jnp i / y \ / IT 

A = — ; sm 1/2^/u/np sm 2y/npu - - 

^cjq 1 . ( 'n'Rg \ . f27cR(j 71 



A = ^ir:^— ^sin ^^^^ sin ^^^^ - - , (13) 
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where we have rewritten equation (11) containing u = K I /2 in terms of the ratio of the 
semi-classical orbital radius Rg and the modulation period a. 

To better understand the increase in amplitude in the magnetically modulated graphene 
system compared to the electrically modulated one we consider the difference in bandwidths 
in the two cases 13|, Il5|. Important feature is the additional sin ( ^^ j term and a factor of 
^/n in the perturbed energy eigenvalues for the magnetically modulated case which is absent 
in the electrically modulated case. For large value of n, sin ( ^^^ j can be taken to be equal 
to ^^^ and the amplitude of the oscillations in the bandwidth of the magnetically modu- 
lated graphene becomes smaller due to this factor compared to the magnetically modulated 
2DEG. The result is that the bandwidth in the magnetically modulated (~^\/^~ 1-5 ) 
case is approximately greater by a factor of 1.5 compared to the electrically modulated 
j — j^= = 0.87 J graphene system at magnetic field strength B=l Tesla and it is smaller 

by a factor of 5 compared to the magnetically modulated (^2^ = 7.6) standard electron gas. 



B. Classical description 

We now give a classical explanation of the asymptotic expression of bandwidth obtained 
in equation (13) which is essentially a large n expression. The classical equations of motion 
along the x and y directions are x{t) = Xq + RgSin(ujgt + ip) and y{t) = yo + RgCosiugt + (/?), 
respectively, where ujg = vv2eB, is the cyclotron frequency for Dirac electrons, Rg is the 
cyclotron orbit radius in graphene while Xq and yo are the center coordinates and (/? is the 
phase factor. Note that Uc = — for standard electrons which is not the same as what we have 
for Dirac electrons given above. Now without loss of generality we may take ip = 0. Thus 
the increase in the average energy of the cyclotron motion due to the magnetic modulation 
is evaluated as 

^ 2 



AE{xo) = — uJosm{Kx{t))dt 
to J 



2 



where to is the period of the orbit. Substituting x{t) yields 

AE{xo) = cuoMKRg) cos(irxo) (14) 



with Jo{z) the Bessel function of zero order. For -^^—^ > l,one can approximate the Bessel 
function Jq by a cosine function as follows 



with the result 



^»i^)-(^)'-^"'' ' 



^E[x^)=uA-—\ sin ^-T , (15) 



Tx'^Rg) \ a 4; ' 

which is the classical expression for the bandwidth. This is the same as obtained in equation 
(13) in the limit of large n which is to be expected as the large n limit corresponds to classical 
results. 

V. CONCLUSIONS 

In this work we have analyzed the electronic spectrum of graphene subjected to a magnetic 
field perpendicular to the graphene layer and a unidirectional periodic magnetic modulation. 
We have determined the density of electronic states and the bandwidth of this system. We 
have also considered the asymptotic expression of bandwidth and its relation to a classical 

9 



description. To highlight the effects of modulation on the density of states and bandwidth, we 
have plotted these quantities as a function of the magnetic field for experimentally relevant 
parameters. We have found that oscillations in the density of states and the bandwidth 
are out of phase and larger in amplitude compared to the electrically modulated graphene. 
We also observe that these oscillations are in phase and smaller in amplitude than those of 
magnetically modulated standard electron gas system. 
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